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Abstract. Using a renormalization approach, we study the asymptotic hmit 
distribution of the maximum value in a set of independent and identically distributed 
random variables raised to a power qn that varies monotonically with the sample size n. 
Under these conditions, a non-standard class of max-stable limit distributions, which 
mirror the classical ones, emerges. Furthermore a transition mechanism between the 
classical and the non-standard limit distributions is brought to light. If q„ grows slower 
than a characteristic function g*, the standard limit distributions are recovered, while 
if Qn behaves asymptotically as Aq*, non-standard limit distributions emerge. 

PACS numbers: 05.10.Cc, 02.50.-r, 05.40.-a 
1. Introduction 

Extreme value statistics, that is the statistics of the largest value in a set of random 
variables, has attracted a lot of attention in many different fields, from probability theory 
PP to physics -where fecund interactions with disordered systems [21 El IH 15]; £is well 
as random walks and interface fluctuations [HI El E] have recently flourished- hydrology 
[9], finance [inilll] or engineering [12]. For independent and identically distributed (iid) 
random variables, asymptotic distributions have been known for long [131 El 1121 H]- 
Depending on the tail of the parent distribution (from which the variables in the set are 
drawn at random), three different distributions emerge. For parent distributions with a 
tail decaying faster than any power law, the limit distribution is the well-known Gumbel 
one (which also found interesting applications outside the field of extreme values [T5]). 
If the parent distribution decays as a power law close to infinity, the so-called Frechet 
distribution is obtained, while if it decays algebraically close to an upper bound, the 
WeibuU distribution is reached. 

In spirit, these results bear some similarities with the Central Limit Theorem, which 
addresses a similar issue for the problem of random sums instead of extreme values. 
Interestingly, it has been shown recently that the limit distribution of sums can be 
modified by raising the summed variables to a power that diverges with the number of 
terms in the sum [TB] . Such a problem is actually motivated by the physics of disordered 
systems, as it can be interpreted as the partition function of the Random Energy Model 
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[T7] , one of the simplest disordered models -which led to recent developments in relation 
to extreme value statistics [SlIllH]. This problem also exhibits interesting connections 
to empirical moment estimation in signal processing and multifractal analysis [T8] . 
It is then natural to wonder whether such a procedure, namely raising the random 
variable to a power increasing with the sample size, could generate some non-standard 
distributions as far as extreme values are concerned. In terms of Random Energy Model, 
this would mean considering the statistics of the maximum value of the Boltzmann 
weights (which add up to the partition function). A related, but perhaps more concrete, 
physical example is the statistics of the largest trapping time in a trap model [TH], in 
which particles are trapped in deep energy wells and can escape only through thermal 
activation. These extreme times are known to play an important role in this context. 
In the limits of low temperature and large number of traps, the statistics of the largest 
trapping time could depend on the way the two limits are taken. 

This issue has been recently addressed in the mathematical literature [20j, following 
the work by Ben Arous and coworkers on the problem of sums, obtaining precise 
results about a transition between the Gumbel attraction domain and the Frechet 
attraction domain for a specific class of distributions. In addition, this problem has 
some connections with the question of the existence of different limit distributions using 
power rescaling procedures [211 [221 [231 [211 [25] • 

In the present contribution, we address the general issue of the limit distribution 
of the maximum value in a set of random variables raised to a power exponent 
diverging with sample size, by generalizing the renormalization group approach recently 
introduced to deal with finite size effects in standard extreme value statistics [2S1 123 [2B] ■ 
For exponents increasing as a power law of the sample size, we find non-standard limit 
distributions, which turn out to be related by an exponential change of variables to 
the standard limit distributions. We clarify this surprising relationship using a simple 
argument based on the behaviour of the rescaling factors. 

2. Problem statement 

Starting from a set of iid variables (VFj)i=i,...,n5 we consider the maximum in the 
set, namely: 

Mf = max{iyi,...iy4. (1) 

Classical extreme value theorems yield asymptotic convergence results for as a 
function of the behaviour of the tail of the probability distribution of W . More precisely, 
there exist two sequences a„ and such that the cumulative distribution of the rescaled 
random variable F„ = (M^— / an converges to the hmit cumulative distribution F^^iy) 
defined by 
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r e-i^+Cv)-'"^ if C > {y> -1/C) 
^dv) = e--^ if C = (2) 

[e-i^Hy)-^^< if^<0 (y<-l/C) 

The case C > corresponds to variables Wi with a distribution decaying as a power law 
at infinity, while the case C < rather corresponds to a power-law decay close to an 
upper bound. Finally, the value C = is obtained for distributions decaying faster than 
any power law (either at infinity or close to an upper bound) [131 dH [T21 [I] . 

We wish to investigate whether a n-dependent transformation of the variables Wi 
may lead to asymptotic distributions different from the present ones. We are especially 
interested in power transformations of the form Ui^n = Wfl" where g„ depends on n, but 
begin by considering the general class of transformations 

U,,n = oj-\W,), ^ = l,...n, (3) 

where Un consists in an increasing bijective function 0. One can express the cumulative 
distribution of ?7 as : 

Fu,niu) = FwiuJniu)). (4) 

In the power transformation case, this definition of u;„(m) leads to 

cj„(m) =M^/''", u>0. (5) 

This transformation is reminiscent of the study in [16j concerning the behaviour of sums 
of random exponentials and notably the failure of the classical central limit theorem for 
rapidly growing powers. 

On the one hand, if Un varies sufficiently slowly as a function of n, it is expected 
that the transformation Un does not affect the limit distribution of the maximum. 
On the other hand, for well-chosen transformations, it should be possible to attain 
new types of limit distributions. Let us consider the transformed maximum = 
max{f/i j^, . . . Un,n}- The cumulative distribution Fnirn) of can be expressed in terms 
of the cumulative function Fw of the variable W as : 

F„(m) = Fvy(u;„(m))". (6) 

In the following section, we devise a renormalization group formulation of Eq. ([6]), which 
allows us to derive in a straightforward way the possible fixed point distributions. 

X Considering decreasing bijective functions is another possibility which would require only minor 
changes. 
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3. Renormalization approach 

3.1. Renormalization transformation and standardization conditions 

Following Refs. [271 [2H], we introduce the functions (7n(^) = ^ ln[— In F„(m)] as well as 
gw{w) = — ln[— InF-iy(w)], and recast Eq. ([6]) into the form 

Qnim) = gw{uJn{m)'j - Inn. (7) 

As in the case of standard convergence theorems, it is useful in order to converge 
to a non-degenerate limit distribution to rescale the maximum value through 
Xn = {M^ — 6„)/a„, where a„ and 6„ are chosen so as to meet some specific conditions 
(for instance, fixing the values of the first two moments). In addition, it is also convenient 
to consider as a real variable rather than an integer one and to define the variable 
s = Inn. We thus assume that a;„(m) can be extended to real values of n, and we define 
the function u{m,s) = uJe!>{m). Altogether, one obtains from Eq. the following 
evolution equation, in terms of the variable X = X^s: 

g{x, s) = gw{uj{a{s)x + b{s), s)) - s, (8) 

where exp[— ex]i{—g{x, s)] is the cumulative distribution of X. In order to determine a{s) 
and b{s), one needs to impose 'standardization' conditions on g{x,s). Such constraints 
are arbitrary to some extent, and may differ depending whether one is interested in 
practical problems or in theoretical approaches. In practical applications, fixing some 
moments of the distribution (e.g., the first two moments) may be convenient. In contrast, 
it turns out that for theoretical purposes, fixing the value of g{x, s) and of its derivative 
dxg{x, s) at a given value of x is an easier condition to implement We thus choose 
the same conditions as in [2H], namely 

9iO,s) = 0, (9) 

dxg{o,s) = 1. 

These conditions imply 

gw{uj{b{s),s)^ = s, (10) 
dsUj{b{s),s) 

a{s) = b{s) + - — , 11 

dmCd{b{.s),s) 

with b'{s) the derivative of b{s) and dm<^ indicates the derivative with respect to the 
first argument of u{m,s). In order to simplify the expression of gw, we contract the 
transformation ijj{x, s) and the rescaling operation into a single transformation T{x, s): 

T{x, s) = u{a{s)x + b{s), s), (12) 

which leads to 

g{x, s) = gwiT{x, s)) - s. (13) 
§ Throughout the paper, we use the notations dx = d/dx and ds = d/ds. 
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3.2. Partial differential equation for the flow 

The functional equation ( |T3l) can be converted into a partial differential equation. We 
first differentiate g{x,s) with respect to x and s: 

d^g{x, s) = g'wiT{x, s)) d^T{x, s), (14) 
dsgix, s) = g'w{T{x, s)) dsT{x, s) - 1. (15) 

Reinjecting Eq. ( fT4l) into Eq. ( |T5i) to eliminate g'^r, we obtain 

dsg{x,s) =U{x,s)d:,g{x,s) -I, (16) 

where we have defined 

The function U{x,s) can be expressed explicitly using Eqs. (IT^ and (IT7|) . We start by 
computing c^^T and dsT (in order to lighten the notations, we drop in the following the 
explicit s dependence of the parameters a{s) and b{s)): 

dxT{x, s) = adm^^{ax + b, s), (18) 
dsT{x, s) = {a'x + b') dm^^{ax + b,s) + dsu{ax + b,s). (19) 

We finally obtain 

a' b' 1 dgijjiax + b, s) , , 

U{x,s) = -x + - + -^^^ (20 

a a aOm^[ax + o,s) 

From now on, we focus on the case of a power-law transformation uj{m,s) = nn}/'^^^\ 
with m > 0. One has 

q{s)m 



dmUj{m, s) = —— — , (21) 
q[s)m 

q'(s)\nm 

OsUj{m,s) = ;r^-r — u{m,s), (22) 



and U{x,s) reads 



a' b' a' ( b\ 
U{x,s) = —x^ — x+- ln(ax + 6). (23) 



Taking into account the standardization condition (fTTi) . which reads, in the case of 
a power-law transformation w(m, s). 



a = h'--b\nb, (24) 



one can rewrite Eq. fl25]) as 



U{x, s) = ( ^ + (9, In ^ ) X + 1 - M X + ^ j In ( ^x + 1 ) . (25) 
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If one defines 

A(^) = ^, 5{s) = dM\{s), ^i^s) = \{s) + 5{s), Q(s) = ^, (26) 
b{s) q{s) 



it is possible to rewrite Eq. ( 1251) in a more compact way as 

U{x, s) = l + ^{s)x - Q{s) {X{s)x + 1) + (27) 

3. 3. Fixed point distributions 

A stationary solution, that is a solution of Eq. ([TB]) satisfying dsg{x, s) = 0, can be 
obtained on condition that U{x, s) is independent of s, namely U{x, s) = lA{x). In this 
case, the stationary solution g{x) is determined by integrating the differential equation 

</i^) = (28) 

with the condition (^(0) = 0. 

We now investigate under which condition U{x, s) becomes independent of s. 

3.4- Case Q = 0; recovering standard limit distributions 
In the case Q = 0, one has 

U{x,s) = l + j{s)x. (29) 
The condition dgU = yields that 7(5) must be equal to a constant 7. Consequently, 

g'ix) = (30) 

1 + 7X 

Taking into account the standardization condition we obtain the fixed point function 

g(x) = -ln(l + jx). (31) 

7 

Reformulating this result in terms of the cumulative distribution F{x) = 
exp[— exp(— 5f(x))], one recovers the standard limit distributions given in Eq. ([2]), 

F(x) = J'^(x) = exp [-(1 + 7x)~^] , 1 + 7X>0, (32) 

where 7 plays the role of the parameter (. Hence classical limit distributions are 
retrieved in the case where the power q{s) is a constant, namely q{s) = go- This 
result was expected: if X belongs to the attraction domain of J-"^, X'^° either belongs to 
the attraction domain of J^^/go for > 0, or to the attraction domain of J-",^ otherwise. 
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3.5. Case Q ^ 0; emergence of non-standard stable distributions 

We now turn to the case Q 0. From Eq. ( 1271) . it is clear that U{x,s) is independent 
of s if Q{s) and A(s) are constants, Q{s) = Q and A(s) = A > 0. Hence stationary 
solutions only exist if q{s) is of the form q{s) = K e^^ , with K > Q a, real constant. 
Inserting Q and A in Eq. ( |25l) . Ll{x) takes the form 

W(a;) = 1 + Aa;-^(1 + Ax)ln(l + Ax). (33) 
A 

Combining Eq. ^ with Eq. ^ leads to 

9\x) = J ^- ^. (34) 

(1 + Aa;) (l-f ln(l + Aa;)) 

By definition, g{x) has to be an increasing function of x so that g'{x) > 0, which implies 
that X belongs to a restricted range of values, Xmin < a; < Xmax- The bounds Xmin and 
a^max are determined by the conditions 1 + 7X > and 1 — y ln(l + Ax) > 0. Assuming 
Q > 0, one finds 

Xmin = -J , W = J (e^/^^ - 1) , A > 0. (35) 
Similarly, for Q < 0, 



Xmin = J (e^/^ - 1) , Xniax = +00 , A > 0. (36) 

Coming back to Eq. fl51|l . this equation can be integrated into 



aix) = --^ In ( 1 - — ln(l + 7x) ) , x^m < x < x^ax, (37) 

also taking into account the condition g{0) = 0. The corresponding cumulative 
distribution J'x,q{x) reads 



J^^Q^x) = exp 



^ln(l + Ax) 
A 



which generalizes the standard extreme value distributions; for Q = 1, this expression 
reduces to a power law on the interval [xmin, a;max] ■ Note that the expression (l38l) of 
J^x,q{x) converges, in the limit Q — > 0, to the standard distribution J-'\{x) given in 
Eq. (1321) . In addition, it is interesting to note that in the limit A — >• 0, the cumulative 
Tx,Q converges to J^_q. 

Remarkably, these generalized extreme value distributions are closely related to the 
standard ones. In term of cumulative, expression (1551) implies that 

J^x,q{x)=T-q{H1 + Xx)/X). (39) 
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Figure 1. (Colour online) Probability density function Px.q{x) associated to the 
generalized extreme distribution T\,q{x). Top left: Q = 0.7, solid black line A = 0.5, 
dashed red line A = 1.2. Top right: (5 = 1, solid black line A = 1, dashed red line 
A — 0.4. Bottom left: Q — 1.5, solid black line A = 2, dashed red line A = 1. Bottom 
right: Q — —1, solid black line A = 1, dashed red line A — 0.01. 



Incidentally, setting ( = —Q, Eq. (1391) implies that if Z is a random variable with 
cumulative J-'<j, then the cumulative of [exp (AZ) — 1]/A is J^\-(- Note that Eq. (!39il 
is also valid for Q = 0, and corresponds in this case to a known relation between the 
Gumbel and the Frechet distributions. 

It should also be noted that for Q > these limits distributions are equivalent, 
up to an affine transformation, to the non-standard-distribution obtained in [2S] using 
power transformations as rescaling factors. This suggests that q{s) plays a role akin to 
a rescaling factor. However, let us emphasize that our approach is different in spirit 
from that of Ref. [25]. In the latter, the power g„ is used as an adjustable rescaling 
parameter allowing the transformed distribution to converge to a non-degenerate limit. 
In the present paper, we consider the power as a given function, and we study the 
non-degenerate limit distributions that are obtained through a standard affine rescaling 
of the data. It is then a priori not obvious that the same stable distributions should 
emerge from the two procedures. 

Deriving the probability density function P\.q{x) of these generalized extreme value 
distributions yields the expression 



PX,Q[X) 



1 + Xx 



^1 n 

-ln(l 



Ax) 



exp 



l-^ln(l + Ax) 
A 



(40) 



As depicted by Fig. ([T]), this family of probability density functions has a non-trivial 
behaviour near its bounds. For Q > 0, the support of pa,q is (— 1/A, [exp{X/Q) — 1]/A). 
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If we consider the asymptotic behaviour in the hmit x 
obtain 



Q 
y 



ln(6) 



^ |Slne|i/ 
gl A 



1/A + e, with e — 0, we 

(41) 



Therefore, we have a crossover at Q = 1 where the probabihty density corresponds to 
a power law of exponent 1/A — 1. For Q > 1, Px,q{x) diverges when x -> —1/A faster 
than any power functions, whereas for Q < 1 and all n, Px,q{x) and all its derivatives 
converge to 0. 

Similarly, in the neighborhood of the upper bound x = [exp(A/Q) — 1]/A — e, we 

have ^ 

-VQ rn^-A/Q^^Q-^ (42) 



P\,Q la;max - ej 



and a singularity appears at the upper bound of the distribution for Q > 1. 

Conversely, for Q < 0, {exp{\/Q) — 1)/A becomes the lower bound of this 
distribution, and using x = Xmm + e and e — > leads to 



P\,Q [Xmin + e) 



-X/Q 



{\Q\ 



g-A/Q \Q- 



exp 



(43) 



Since Q < 0, the exponential term is dominant and the probability density function and 
all its derivatives converge to in the neighborhood of x^^i^. On the other hand, for 
X — > +00, F{x) behaves as the exponential of the power of a logarithm: 



F{x) 



x— >-+oo 



exp 





— Inx 


1 ' 

Q 




A 





(44) 



There is therefore four distinct asymptotic behaviours for the function px^q{x) depending 
on the value of Q {Q < 0, < Q < 1, Q = 1, Q > 1). The function px q{x) is shown in 
Fig. [T] for several values corresponding to these distinct domains. 



4. Attraction domain of the non-standard limit laws 



Having found non-standard asymptotic forms of the extreme value distribution in 
Section [31 we now consider their domain of attraction and we use the partial differential 
equation of the flow to develop a heuristic description of these attraction domains. This 
description provides a good approximation of the standard attraction domains and gives 
us precise insights into the non-standard attraction domains. Finally, these heuristic 
results are confirmed in section fusing an independent approach which sheds some light 
on the relationship between standard and non-standard laws. 

Although the stationary equation resulting from (1251) and (1551) cannot be used 
directly to describe the dynamic of the general equation ([TB]) . it underlines at least the 
existence of two important control parameters in the stationary case, namely Q and A. 

If we return to the general flow equation (1161) and the general expression of U 
described by Eq. ( 12 7p . these two variables can be interpreted as an external forcing. 
Heuristically, the convergence of Eq. (fT6l) should be driven by the asymptotic behaviour 
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of these parameters. Notably, if 7(s) and Q{s) converge simultaneously toward the 
respective finite limits F and Q, it is plausible to expect the transformed maximum to 
converge toward the limit law Jt,q- As a corollary, for laws belonging to a classical 
domain with parameter ( in the case Q{s) = 0, one should obtain F = which sets an 
interesting test case for our heuristic argument. 



4-1- Asymptotic behaviours of the forcing parameter 

If the term Q{s) clearly denotes the external forcing due to the power transformation, 
the interpretation of A(s) and 5{s) is hazier. From its definition. 



A(s) = ds\nb{s) ~q'{s)lnb{s) 
^ 9(g) 



(45) 



It is thus possible to factorize A(s) into the factor q{s) and a term Ao(s) depending 
only on the cumulative function of the parent distribution. Using the change of variable 
s = g{x) and the relation g{x) = — ln(— InF(x)) leads to 

\F{x) lnF{x)\ 

M9{x)) = ^^.(^) • (46) 

One should note that Ao(s) > due to the properties of the cumulative function. 
Moreover, defining the complementary cumulative function F(x) = 1 — F{x) leads to 
the asymptotic expression 

Xo{g{x)) ~ (47) 

The variable Ao(s) corresponds to the inverse of the local power exponent of the 
complementary cumulative function F at point x = g{s). If Ao(s) admits a non-zero 
finite limit Aq at +oo, there exists a normalized slowly varying function such that 

F{x) = L{x)x'^/^". (48) 

In other words, the parameter Aq describes the power behaviour of the tail of the 
distribution. Since for Q{s) = 0, A(s) and Ao(s) are identical, laws belonging to 
the classical domains of attraction provide enlightening examples of the asymptotic 
behaviour of Ao(s). Notably, if a law belongs to the Frechet attraction domain of 
parameter (, there is a slowly varying function L such that 

F{x) = L{x)x-'~ . (49) 

II The slowly varying functions class gathers constants, logarithms and all functions satisfying the 
asymptotic relation: Va > 0, L[ax) ~ L[x). A slowly varying function is said to be normalized 

if lima;_>+oo xL'{x)/L{x) = 25J p. 15]. Since we are only interested in the generic properties of slowly 
varying functions, we most often use the same notation L for any function belonging to this class. 
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With the regularity assumption that L is normahzed, we obtain that Aq = C- Conversely, 
for a law belonging to the WeibuU domain, we have a slowly varying function L such 
that 

F{x) = {xf - x)-^I^l({xf - x)"^) (50) 

where xp is the end point of the distribution. Under the assumption that L is normalized, 
we have 

The Gumbel domain of attraction is harder to characterize. However, the "exponential 
power" laws 

= exp(-x") (52) 

often constitute an interesting subset of this domain. For such laws, a short calculation 
leads to 

\,{g{x)) ^ ^. (53) 

In the general law belongs to the Gumbel attraction domain, if and only if 

1/F{x) is a rapidly varying function. Using an integral characterization of rapidly 
varying functions [2^ p. 178], it is possible to show, that under a regularity assumption, 
for any law belonging to the Gumbel attraction domain, Ao(s) — )■ 0. Consequently, for 
laws in the Weibull and Gumbel attraction domains, we have Aq = 0. In summary, for 
laws belonging to a standard attraction domain of parameter (, one has Aq = max(0, (). 
Similarly, 6{s) can be decomposed into 

S{s)=Q{s)+So{s), do{s) = dJnXois). (54) 

One interesting consequence of the previous equation is that the asymptotic behaviours 
of Sq{s) and Ao(s) are entwined. If we assume that both Sq{s) and Xo{s) admit a finite 
limit when s goes to +oo, with lims_j.oo Ao(s) > 0, then one necessarily has lims_j.oo 6o = 0. 
This situation corresponds to the Frechet class. Conversely, if lims_s.oo Ao(s) = (a 
typical situation in the Gumbel and Weibull classes), and if So{s) admits a finite limit 
Aq when s — j- cxd, then Aq < 0. After some tedious calculations, one can expand the 
expression (IMl) of So{s) as 

So{g{x)) = I lnF(a;)| — 1 - Ao(^(a;)) 

(55) 



Going on with our study of the classical domains, we consider a law belonging to the 
Weibull attraction domain with parameter ( and slowly varying function L. Under the 
regularity assumption that both L and L' are normalized slowly varying function, it is 
possible to verify that Aq = C- The same computation for the "exponential power" laws 
(a typical example of laws belonging to the Gumbel domain) leads to 

5oigix)) ~ (56) 

x— )-+oo 
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and therefore, Aq = 0. In other words, the hmits Aq and Aq are respectively the 
positive and negative parts of Tq = hms^oo('^o('5) + Ao(s)), namely Aq = max(0, Tq) and 
Aq = min(0,ro). For a law belonging to the standard attraction domain of parameter 
(, we have as expected Fq = C- So, for Q = 0, we have recovered the classical results 
concerning the attraction domain (ignoring the difference in regularity assumptions). 
And from this point, it is easy to extend this convergence result to the non-standard 
case (Q 7^ 0). 

Piecing together Eq. (HSi) and 7(5) can also be rewritten as 

^{s) = 6o{s) + Q{s)+q{s)\o{s). (57) 

The previous expression outlines an interesting interplay between the tail behaviour 
of the parent distribution and the power transformation. The tail behaviour of the 
distribution is responsible for the parameter 6q whereas Q{s) is directly derived from 
the choice of q{s). The term q{s)Xo{s) represents the interaction between the two effects. 



4-2. Changing the attraction domain by varying q{s) 

If we consider a fixed parent law, the functions Sq{s) and Ao(s) are then fixed. Therefore 
the only free parameter is q{s). In this situation, the limit F is determined by the 
asymptotic behaviour of A(s) = g(s)Ao(s). Let us introduce the limit A = lim5^ooA(s), 
when it exists. If we define a characteristic power scale by 

q*is) = (58) 

we obtain 

A = lim (59) 
q*[sj 

In other words, taking the limit s — ?■ 00 in Eq. (1571) . one finds that F is only dependent 
of Aq, Q and of the limit ratio between q{s) and the characteristic power scale q*{s). 
From this point, two different situations arise. On the one hand, a special case appears 
if lim^ g(s)/g*(s) = 0. Indeed, if lims_j,oo A(s) = 0, the expression ofU{x,s) can be 
linearised as 

U{x,s) ^l + -f{s)x -Q{s)x. (60) 

So if 7(s) and Q{s) admit the finite limits F and Q, we obtain a partial differential 
equation corresponding to a standard limit distribution with parameter ( = Aq. This 
implies that the WeibuU and Gumbel domains are unaffected by such a transformation. 
However, for a law belonging to the Frechet domain of parameter (, we have Aq = C ^■iicl 
A = 0. This means that A = is only possible with limg q{s) = and Q < 0. In other 
words, a decreasing power destabilizes all the Frechet domains and any law belonging 
to a Frechet domain will converge towards a Gumbel distribution once exposed to a 
vanishing power transformation. 
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On the other hand, non-standard hmit laws appears when q{s) ~ Aq*{s). In this 
case, we have 

r = Ao + Q + A. (61) 
Moreover, the relation Q{s) = q'{s)/q{s) leads to 

Qis) = -5ois) + (62) 
A(s) 

So, if we assume that lims_j,_|_oo A'(s)/A(s) = 0, one obtains that 

(63) 

r = A. 

So using a power transformation with q{s) ~ Ag*(s) leads to the non-standard limit 
law J-'x-Ao- Moreover, it is possible to compute an analytic representation of q* by 
defining the error term e(s) = Q{s) + Aq. Using Eq. 0261) leads to the exact differential 
equation 

ds\nq*{s) = -Ao + e{s). (64) 

Solving this equation results in 

q: = L*{n)n-''\ (65) 



L*(n)=exp(/;iiMdt^ 



(66) 



The factor L*{n) corresponds to a corrective term depending on the fine convergence 
structure of Sq{s) and A(s). Moreover, since e(lnn) — > 0, Eq. fIBB]) corresponds exactly to 
the integral representation of a normalised slowly varying function. In other words, L*{n) 
is a slowly varying correction to the power behaviour of g*, if Aq 7^ 0. However, if Aq = 0, 
this slowly varying term becomes preponderant. Notably for Frechet distributions, 
Ao = C implies that g* ~ On the other hand, for the Gumbel attraction domain, 
Ao(s) — )■ implies that L*{n) — > +00. For instance, in the case of the exponential power 
laws, we have So{s) = e(s) ~ ^, therefore g„ Inn {K > 0). Computing Ao(s) yields 
the more precise expression 

g* ~ a Inn. (67) 

Consequently, using the right power transformation (g„ = AL*(n)n''''), a law 
belonging to the Weibull attraction domain can be forced to converge into a non- 
standard limit distribution J-a,-Ao- On the other hand, after applying a slowly varying 
power transformation (g„ = Aa In n for the exponential power laws) to a law belonging 
to the Gumbel attraction domain, one can force the convergence towards a Frechet 
limit law. This is perfectly consistent with the (more precise) result obtained in [20] 
concerning the transition between the Gumbel domain and the Frechet domain for a 
specific power transformation. 

However, considering parent distributions from standard domains only leads to 
non-standard distributions with positive Q. The domains of attraction of non-standard 
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JTclIcllL idW 


Functional dependence of g„ 


L{n)n9, Q < 


A 


Ld{n) 


L{n)n^, Q > 


VVt;lL)Ull LlUlli. 

C<0, L*{n) 


_F, 0=ICU = AL- ^^^^ 


vJU-lliUt-l U-Ulli. 

L*{n) 




L = XL* 


None 


Frechet dom. 

C>o 






None 


Log. powers 
a > 0, L*{n) 


Q = -a, L = XL* 
J^a None 



Table 1. Classification of the limit distributions according to the functional 
dependence of the exponent g„. The function Ld{n) is a slowly varying function that 
goes to infinity with n. The slowly varying function L*{n) characterizes the parent law 
-see Eq. (|65p . 'None' means that no (non-degenerate) limit distribution emerges. 



Parent law 


Limit of the ratio qnl<in 


Qn/q*n 


Qn/ql ^ A > 




Frechet domain, C > 


Gumbel J^q 


Frechet J-a 


None 


Gumbel domain, J^q 


Gumbel J-'q 


Frechet J-a 


None 


Weibull domain, C < 


Weibull J"^ 


Non-standard ^a,|c| 


None 


Logarithmic power, a > 


Frechet Fa 


Non-standard J-'a^^o 


None 



Table 2. Classification of the limit distributions according to the limit of the ratio 
Qn/Qn when n ^ oo. 

laws with negative Q correspond to parent laws which do not belong to the standard 
attraction domains but can be "renormalised" using the decreasing power law. One 
example of such laws would be the logarithm-power law with 

= 1 -L(lna;)(lna;)--, (68) 

where both L and its derivative L' are normalized slowly varying functions and a > 0. 
A short calculation shows that XQ{g{x)) ~ alnx and Aq = a. Thus these laws belong 
to the attraction domain of J-'x-a with a > and A > 0. The associated power scale 
correspond to a decreasing power functions g* = L*(n)n^". Going further, with a power 
transformation decreasing faster than these logarithm-power law converge towards 
a Frechet distribution of parameter a. 

In conclusion, our heuristic analysis brought to light the existence of transitions 
between the standard and non-standard attraction domains when modifying the n- 
dependence of the power transformation g„ (or q{s)). As a byproduct, we also obtained 
that power transformations can be used to "renormalize" laws beyond the standard 
attraction domains (a typical example being logarithmic-power laws) in such a way that 
they converge to a non-degenerate distribution. This set of transitions is summarized 
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in Table [TJ A somehow more compact presentation can be obtained by using the hmit 
of the ratio Qn/Qn order to classify the limit distributions, as seen in Table [21 

5. Alternative approach to the characterization of non-standard attraction 
domains 

5.1. Rescaling factors for transformed maximum 

In the previous section, heuristic arguments have shown that only very specific 
transitions are possible between the standard and non-standard attraction domains. In 
order to verify this result, this section will present an alternative approach to the study 
of the transformed maximum in which we shall tie back the convergence behaviour of 
the transformed and the non-transformed variables. 

Let us consider a random variable W which belongs to the domain of attraction of 
the limit law J-"^. This means that there exists a renormalization sequence such 
that 

Vx, Fw{anX + PnT ^ ^d^)- (69) 

These renormalization sequences are well-known and expressions are available in the 
literature [TIUfT]. 

Similarly, the w-transformed variable converges in distribution if and only if 
there exist F and a sequence (a„, 6„) such that 

Vx, Fy/{uJn{anX + hn)] ^HdcJ^ix). (70) 

However, in this extension of the previous problem, no general conditions of 
convergence are known and the choice of the renormalization sequence becomes more 
difficult. We propose in the next section to exploit the striking similarity between the 
two previous equations to obtain convergence conditions for the power transformation 
in specific cases. 

First, one can remark that the only difference between the two previous equations 
lies in the term (a„x -|- /3„) in Eq. (jSH]) which becomes aj„(a„x + 6„) in Eq. (1701) . 
This similarity suggests a simple way to obtain the convergence in distribution of the 
transformed maximum by exploiting our knowledge of the renormalization factor of the 
original distribution. If there exists a renormalization sequence (a„, fen) such that 

VX, Un{anX + hn) ~ ttnl'^x) + (5n, (71) 

n— >-+oo 

and further assuming that 

lim Fw{oo{anX + hn)T = lim Fw{any{x) + l3nT , (72) 



one obtains 

jr(x) = 7q{v{x)). (73) 
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Therefore, if the conditions (172]) and (ITTIl are satisfied it is possible to link the 
transformed limit distribution and the standard limit distribution. 

Condition ( 1721) corresponds to a quite technical convergence problem. For now, we 



assume that this condition is satisfied. We show in Appendix A that this condition holds 
for our proposed choice of (a„, Condition (17T1) is more interesting and through z/(x) 
defines the kind of transition. For a power transformation, it can be read as 

Vx, {anX + bnY^"" ~ a„z/(x)+/3„. (74) 

n— f +00 

Considering Table [2], four distinct transitions should be possible. For laws belonging 
to the Gumbel or Weibull domain, one should have either i'{x) = x for lim„_!.oo Qn/Qn — 
or v{x) = ln(l + Ax) /A for limn^oo Qn/Qn = A > 0. Similarly, in the Frechet domain 
with parameter (, transitions should appear for exponents g„ converging to a finite 
value, considering u(x) = [(1 +AxY^^ — !]/( (with A = lim„_j.oo qn/C)-, and for vanishing 
exponents qn considering z/(x) = [exp(Ca;) — 1]/C . 

So, using the insight gained from section HI it is natural to study separately the 
behaviour of (a^x + fen)^''''" for diverging, vanishing and converging power 

5. 2. Diverging powers qn 

In the case where hm„_j.oo <?« = +oo, the expression (a„x + &„)^/''" has two distinct 
asymptotic behaviours which lead to two different convergence regimes. 
First, if ttn/hn — > 0, one has the asymptotic behaviour 

n— >+oo 

{anx + hnf''- ~ h]!'^- (l + -^x") . (75) 
The condition Eq. (I74|) is then satisfied if 

Pn 

_ , On (76) 

n "nqn p. 

Pn 

This choice of (a„, 6„) is compatible with the assumption a„/6„ — if 

n— s>+oo 

Tt 0- (77) 

The term an/ Pn corresponds to the parameter \o{n) defined in section jH So using 
the definition of the characteristic exponent g* = l/Ao(^) introduced in Eq. (15S|) . the 
previous result states that if g„ is negligible with respect to g* (that is, lim„_j.oo q-a/qn = 0) 
then the transformed maximum converges towards the same limit distribution as the 
original maximum. This is the expected result from section |H 

The second asymptotic behaviour arises when a„/6„ is a constant. By factorizing 
one obtains 
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(a„x + 6„y/'?" ^ 5^-^. 1 + ^ ^ . (78) 



n—>-+oo 



Eq. dH]) is satisfied if: 



^ -> A 



(79) 



/7* n— ^+oo 
- 

In tliis case, z/(x) = (1 + Alnx)/A implies tliat X„ converges towards tlie non-standard 
limit laws. More precisely if g„ is asymptotically equivalent to Ag* then the maximum 
converges in distribution towards ^a,-c- Once again, we recover the results of section H] 
for the WeibuU and Gumbel domain. 

5. 3. Vanishing powers g„ 

The next interesting transition appears for vanishing moments. From Table El we know 
that an exponential term should appears in (a„x + bnY^'^"- The easiest way to obtain 
this term is to assume that = hnQnC- Then we have 

~ (6„)^e^^. 

n— >+oo 

Substituting z/(x) = [exp{(x) — 1]/C) one may also satisfy Eq. ( 1711) by assuming 

bn = iPnr 

< «n ^ 1 (80) 

It is therefore possible to go from the Frechet domain to the Gumbel domain using 
any decreasing power transformation, as expected from our heuristic analysis in the 
previous section. This confirms that the Frechet domains are very unstable under power 
transformation. Any vanishing power is enough to change a distribution belonging to 
the Frechet domain to converge towards the Gumbel distribution. 

5.4- Converging powers Qn 

In the case of converging powers, one expects to observe only transitions between Frechet 
domains with distinct parameters C,. Considering the possible translation and dilation, 
one can define without loss of generality = A6„ and 6„ = Then, one obtains 

(a„x + M''''" = /3n(l + Ax)i/^", (81) 
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and choosing z/(x) = [(1 + Ax)'^/^ — leads to the condition 

(82) 




We observe as expected a transition between the Frechet domain of parameter ( and the 
Frechet domain of parameter A. If the choice of z/(a;) can appears to be quite arbitrary, 
it should be noted that it is merely a consequence of the representation chosen for the 
Frechet limit laws. A different choice of representation {F{x) = exp(— leads to 
the far simpler iy(x) = x^^^ . However, it is compelling that our renormalization methods 
have allowed us to shed light on this transition even in this convoluted settings. 

With this transition between Frechet domains, we have recovered all the possible 
transitions from standard attraction domains to other domains described in Table |2] 
using only the insight obtained from our analysis of the partial equation of the 
renormalization flow and elementary arguments on the renormalization coefficients. As 



shown in Appendix A these arguments lead to a rigorous proof of the convergence of 



the transformed maximum. 



6. Conclusion 



In this contribution, the renormalization approach of the problem of maximum 
developed in [2SII2ZII2B] has been extended to the case where the underlying variables W, 
are subjected to a transformation which depends on the sample size. The reduction 
of the problem of maximum to a partial differential equation turns out to be a rather 
straightforward generalization of the standard case and leads in the case of the power 
transformation f/j „ = VT/", to a quite short categorization of the limit distributions. 

Using this categorization, non-standard max-stable laws mirroring the standard 
limit laws have been brought to light. These new limit laws are closely related to the 
standard ones. However the behaviour of the partial differential equation describing 
the evolution of the distribution of the maximum is more complex and involves some 
intriguing interactions between the rate of growth of the power transformation and the 
tail of the distribution. 

These interactions received further investigations by studying the asymptotic 
behaviour of the forcing parameters appearing in the partial differential equation of 
the flow. Using a heuristic argument, it was possible to recover a slight approximation 
of the standard attraction domain. Moreover, the same argument leads to an interesting 
description of the attraction domain of a non-standard law, illustrating the existence 
of specific transitions between the classical limit laws and their mirrors laws, when 
varying the functional dependence of the power g„. These transitions are associated 
with a characteristic power scale g*. If q-n/Qn ~^ power transformation g„ is 

too slow to infiuence the convergence of the maximum towards the standard limit 
distributions. In contrast, if g„ ~ Ag* the distribution converges towards a non-standard 
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limit distribution. Using insights gained from the partial differential equation of flow, it 
was then possible to confirm the existence of these transitions, and to investigate their 
mechanisms using a more direct approach based on a study of the rescaling factors. 
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Appendix A. Convergence problems 

The convergence results obtained in section O are dependent on the condition given in 
Eq. (ESI). This condition is equivalent to F((a„a; + 6„)i/5")"-F(a„z/(x) 0. This 
is immediately true if gw is uniformly continuous. However, this condition is superfluous. 
Let us define 

Xn = anU{x) + f3n, (A.l) 

an 

It is then possible to show that for the four cases presented in sectional lim„_j.oo en = 0. 
For diverging g„ and qn/ln ~^ 0> -^l- (EH]) leads to 

en - (--iV^V (^-3) 
In the case g„ ~ \q*, we have 

.„.Ml...)(i-|)+o(^), (A.4) 

Similarly for vanishing g„, Eq. ( ISUl) yields 

en= (g:-i)(e^^-l) + 0(gn). (A.5) 

Since g„ — in this case and Qn l/C? is therefore a vanishing quantity. And lastly, 
for converging moment and g„ ~ Ag* 

6n % ^((l + Ax)^/^"-(l + Axp) + ^-g: (A.6) 

Using g* — 1/^, this confirms that lim„_j.ooe„ = 0. 

Moreover, by construction lim„^oo ""--^(^^n) ^ IR- It is possible to show that 
F{{anX + 6„)^/^")" - F(a„z/(a;) + /^„)" is equivalent to 

hm ^(^:^+ ^"^") = 1. (A.7) 

If Eq. (IA.7P is satisfied, then the convergence in distribution of is ensured. In 
sections Appendix A.l, Appendix A. 2 and Appendix A. 3, we verify that this condition 
holds for all the standard domains. Consequently, the transition described earlier is 
always valid. 
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Appendix A. 1. Gumhel domain 

The Gumbel attraction domain is the harder to characterize. In order to prove the 
convergence in the general settings, we will temporarily use the standard renormalisation 
factors 

f5n = F-\l/n), (A.8) 

an = E{f3n), (A.9) 

where the function E(x), defined as E(x) = i_p(^^^ Jx°° 1 ~ F(t)dt, satisfies: 

lim E'(x) = 

x— >+oo 

^ , l-F{x + rE{x)) (A-10) 
Vr > 0, lim ^ — ^, , = e 

x^+co 1 — F[X) 

One useful property of E is that for any positive real r 

EiE±lE(-)l ^ 1. (A.11) 

E{X) X^ + OD 

Combining Eqs. flA.lOD . flA.llI) and flA.llI) leads to 



E{Xn) = E{anX + I3n) 

= E{(3n+xE{(5n))r^E{(3n), 



(A. 12) 



so that E{xn) ~ ctn- Consequently, a„en/-E(a;„) — )• 0, and using this result with 
Eq. (lAlOll yields 



) = F{Xn H —, rE{Xn)Xn) 



~ F{Xn) 

So Eq. ( 1A.7[1 holds for the Gumbel attraction domain. 



(A.13) 



Appendix A. 2. Weihull domain 

As stated in Eq. ( l50l) . we have for the Weibull domain 

= l[{xf - x)-^) {xf - x)-^/^. (A.14) 
Moreover, using the definition of the renormalisation factor leads to 

Xf - (3n {xf - f3n)F{f3n) 



an «^+°^ F (/3„) 
Hence, with the assumption that L'/L — > we have 



(A.15) 



lim ^^J^ = ^. (A. 16) 
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Therefore, it is possible to show that 

xp - {xn + a„e„) ~ a„x ^1 + . (A. 17) 

The properties of slowly varying functions yield, for n — oo 

F{xn + a„en) = L ( ^ ) {xp - x„ - a„e„)^ 

\xf - Xn - Q;„e„/ 



Hence Eq. flA.7p holds for the WeibuU domain. 

Appendix A. 3. Frechet domain 

Within the Frechet domain, the proof is immediate using an^n/xn which is directly 
implied by lim„_j.oo e = and lim„^oo Q'n > 0- We have from Eq. 



~ x-^Kl(x„) = F(x„). 



As a result, Eq. (1A.7|) is also satisfied for the Frechet class. 
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